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Abstract 

The finite XXZ spin chain with boundaries is studied. We derive the transfer 



o 

[-««. . matrix from the (/-difference equation discovered by Cherednik and construct its 

eigenstates by the vertex operator approach. We point out that the eigenstates 
with no magnetic fields have a symmetry called the turning symmetry. Making 
use of this symmetry we calculate the spontaneous magnetization in the thermo- 
dynamic limit, which is roughly twice as large as that in the half-infinite XXZ spin 



O chain. 
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1 Introduction 

X, 

The XXZ spin chain is the fundamental model in the study of the integrable systems. It 
satisfies the Yang-Baxter equation, which enables us to clarify exactly the structure of 
the energy levels. Further, the Yang-Baxter equation yields an realization of quantum 
group and nowadays its representation theory promotes the analysis of various integrable 
models . 

Encouraged by the success of the conformal field theory ||, |!|, much attention has 
been focused also to the massive integrable field theories. A remarkable result is the 
discovery of the g-difference equation similar to the Knizhnik-Zamolodchikov equation 
H, which is satisfied by the correlation functions in the integrable models. This equation 
is called the quantum Knizhnik-Zamolodchikov (q-KZ) equation and has been studied 
extensively in both mathematics and physics. Solutions of the q-KZ equation are for- 
mulated with various methods, for examples, the Jackson integral ]7[ |8j, the algebraic 
Bethe ansatz |J and the vertex operators (VOs) in the quantum affine algebra [[I], |3|, |TC|j . 
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■'■e-mail: wadati@monet .phys . s .u-tokyo .ac.jp 
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For instance, the correlation functions in the infinite XXZ spin chain are constructed by 
the VOs and satisfy the q-KZ equation. It is also known that the Smirnov's form factor 



is one of solutions of the q-KZ equation R ITT 



The above discussions can be extended to the half-infinite chain with a special bound- 
ary: there exist a boundary condition that preserves the integrability [T2]] . Such a bound- 
ary condition is described by the boundary K-matrix satisfying the reflection equation 
instead of the Yang-Baxter equation. Using this K-matrix and noting that the VOs 
physically correspond to the half-infinite chain M, we can formulate correlation func- 



tions in the half- infinite XXZ spin chain in terms of the VOs ||13|| . Moreover, they satisfy 



the g-difference equation similar to the q-KZ equation [|TJ]]. In the paper this equation is 
called the boundary q-KZ equation. By manipulation of the boundary q-KZ equation, in 
the half-infinite XXZ spin chain, the two-point function is exactly calculated |15j and by 
the explicit representations of the VOs the correlation functions are found in the integral 
expressions 
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Our purpose in this paper is to calculate the correlation functions of the finite XXZ 
spin chain with boundaries by the vertex operator approach. It is obviously impossible 
to analyze the finite chain in the same manner as above, because the VO implies the 
half- infinite chain. However it has been believed that, in the thermodynamic limit, the 
bounded chain is substituted for the half-infinite chain, because the vertex operator 
approach gives the same result as that obtained by analysis of the Bethe equation for 
a finite chain QTJ, [| [16| . Can it be true? We bring forward a counterargument to this 
expectation by noting special bounded spin chains such that they are invariant under 
overturning (called the turning symmetry). The half- infinite chain never satisfies such 
geometric symmetry. The turning symmetry may strongly influence physical behaviors 
of the model. 

To treat the finite XXZ spin chain we present a new idea. We pay attention to the 
resemblance between the boost operator in the boundary q-KZ equation and the transfer 
matrix of the bounded model. In section |2] we introduce the boundary q-KZ equation and 
get the transfer matrix from this equation. Then the eigenstates of the model is made 
of the VOs. In section |3] we introduce the idea of the turning symmetry and compute 
the vacuum energy making use of this symmetry. We also construct the excited states 
using the VOs of another type and compute these energy levels. It is found that the 
structure of energies is similar to the quasi-particle structure pointed out by Faddeev and 
Takhtajan in the XXX model jl7| . In section f| the thermodynamic limit is considered. 
The vacuum is then expressed by the infinite product of VOs. We show that such infinite 
product is interpreted as a huge group of R-matrices and is simplified by the turning 
symmetry. Its asymptotic form is obtained by applying the Baxter's formula on the 
corner transfer matrix ||18|| . In section |^ the spontaneous magnetization at a boundary 
is calculated by the free field representations of the VOs. We find that it is larger than 
that in the half-infinite XXZ chain. 
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2 Boundary q-KZ Equation and Transfer Matrix 



In this section we introduce the g-difference equation called the boundary q-KZ equation 
and solve it by the vertex operator approach. From the boundary q-KZ equation the 
transfer matrix in the finite XXZ spin chain is squeezed out. Then the solution of the 
boundary q-KZ equation corresponds to the eigenstate of the model, which is expressed 
in terms of the VOs. 

2.1 Boundary q-KZ equation 

As an extension of the q-KZ equation Cherednik proposed the following difference equa- 
tion fnji , 

F(Ci, • • • , q-%, ■ ■ ■ , Cat) = ■ ■ ■ , Cw)*XCi, • • • , Cn), (2.1) 

where Tj(d, . . . , ( N ) is the boost operator given by 

VjU: vv) = R j , j - 1 {( j /q\ j - l )---R j>1 (C j /q%)K j (( j ) 

xK^QRjAQ/Of) ■ • • (2-2) 

The R-matrix Ri,j{() is the intertwiner of U q {sl2) for a map Endc(V^ ®Vj). Kj(C) 
and Kj(C) are the boundary K-matrices defined by a map Endc(Vj) fl2fl . They are 
parameterized by r and f such that Kj(() = Kj((; r) and Kj{Q = Kj((; f), respectively. 
Let v e (e = ±) denote the natural basis of V — C 2 . The matrix elements of R{() and 
K((; r) are defined by the relations, 

R(().(v ei ®v e2 ) = J2 R tk{0{^®v e Q, (2.3) 

e[,e' 2 =± 

K((;r).v e = ^Kf(C;rK, (2.4) 

e'=± 

The solution F((x, . . . , (n) forms an iV-fold tensor space V\ ® • • • ® Vjv- It was shown 
1~5| that the solution of equation (|2.1|) corresponds to the correlation function in the 



half-infinite XXZ spin chain with a boundary. In the paper we call the equation (|2^ 
the boundary q-KZ equation. 

The half-infinite XXZ model is formulated in terms of the VOs of the quantum 
affine algebra U q {sl> i ) 0' 01' ^ so ^ ve the boundary q-KZ equation we need the VO 
$(C) = E e =± $ e(C)-^, the dual VO $*(C) related with $*(C) = $_ £ (-g _1 C) and some 
states \W), (W*\ defined by the following relations, 

XX(C)MC)IW) = MC'W), (2.5) 

e' 
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Then the solution of the boundary q-KZ equation is written as 

F(Ci,...,Gv) = (W*|$(Ci)®---®$(Civ)|W0 

= (W\*«(Ci)'''**AOr)\W)(v ei ®.''®v eir ). (2.7) 

ei,...,ejv 

The proof is readily done. Noting that the commutation relation of the VOs, 

^^g(Cl/C2)^(Cl)^(C 2 ) = $e 2 (C 2 )$e 1 (Cl), (2.8) 



and the definition of the states \W), (W*\, one can check that F(Ci • • • , Cat) satisfies the 
following relations, 

Rj,j+i(Cj/Cj+i)F{Ci, ■ ■ ■ j Cat) = F(d, . . . , Q+i, (j, • • • , Civ), (2.9) 
K n (( n )F(( u ...Xn) = ^(Ci,...,^ 1 ), (2.10) 
^(Ci^Cr 1 ,...,^) = F(g- 2 d,...,Civ), (2.11) 

where the order of the spectral parameters corresponds to that of iV-fold vector spaces. 
These relations show that 

t;.(Ci,...,CaWCi,...,Gv) 

= n hJ :K J :<,\ J ,) • ■■R j>1 (( j /q 2 ( 1 )K j (( j ) 

xR l,j((l(j) ■ ■ ■ ■ ■ ■ R N,j((NCj) 

xK j((j)Rj,N((j/(N) ■ ■■Rj,j+l((j/(j+l) 
xF(Ci,...,CaO 
= II,, .(Q/q'Q ,) • ■■R j>1 (Q/q%)K j (( j ) 

x Ri,j(d(j) • • • • • • R N,j((N(j) 

xF{Cu ■ • • > 0-i) 0+i) • • • > Cat, C7 1 ) 

= r m l(0/<? 2 ,) • ■■R J AQ/q 2 Ci)K 3 (Q 

x F((i ) Ci> • • • > O-i) O+i' • • • 5 Cat) 
= /•'((, 1 m/ %,(,,., (,v). (2.12) 

which completes the proof. The boundary q-KZ equation is useful in calculating the 



two-point correlation functions not only in the XXZ model but in the XYZ model [15 



2.2 Transfer Matrix in the Finite XXZ Spin Chain 

One may notice that the boost operator fl2.2|) is similar to the transfer matrix of the 
XXZ model. Now we consider the case such that the parameter q~ 2 does not appear, 

c,v) = ^_ 1 (o/o-i)---^,i(o/Ci)^(o) 

x-Rij(CiO) ' ' ' Rj-ij{Cj-iCj)Rj+ij((j+i(j) ■ ■ ■ R-N,j((N(j) 

xK^QRjMQ/Cn) ■ ■■R j , j+1 (Q/C j+1 ). (2.13) 
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This operator is obtained when the defining relation (|2.6| ) is modified (see (|2.31|) ). It 
is expected that the solution of the boundary q-KZ equation for this operator gives the 
eigenstate of the XXZ model. First we show that the operator fl2.13| ) can be regarded as 
the transfer matrix for the finite XXZ spin chain. We write the matrix elements of the 
R-matrix as follows, 



ittt(C) = JCl(C)=7(0, (2 - 14) 

:i - c 2 )< ? 

1 - g2 C 2 

:i-g 2 )c , 

1 - q 2 ( 2 



RXZ(C) = R--X(C) = V-^7(C), (2-15) 
R- + + _(0 = R±+(0 = T— 07(0, (2.16) 



with the normalization factor 7(C) 

7(0 = 7 r^AvrQ4r » (*;p)oo = f[(i ( 2 -^) 



n=0 



Other elements are 0. The normalization factor ^(Q is derived from the commutation 
relation of the VOs ( |2.8| ). The R-matrix satisfies the Yang-Baxter equation, the unitarity 
relation and the crossing symmetry, 

^12(Cl/C 2 )i?13(Cl/C3)i?23(C 2 /C3) = i?23(C 2 /C3) J Rl3(Cl/C3)i?12(Cl/C 2 ), (2.18) 

J Ri 2 (Ci/C 2 )^ 2 i(C 2 /Ci) = i, (2.19) 

^1i(Ci/c 2 ) = ^::;(c 2 Mi). (2.20) 

The K-matrix is given by the following diagonal matrix, 

K+(C, r) = r), Kz(C, r) = «(C; r), (2.21) 

C -r 

u v (g 4 < 2 ; g 4 )oo(g 2 rC- 2 ; g 4 Ug 6 C 4 ; g 8 )oo(g 8 C 4 ; g 8 )oo (??2) 

K[( ^ r) (g 4 rC- 2 ;g 4 )oo(g 2 rC 2 ;g 4 )oo(g 6 C- 4 ;g 8 )oo(g 8 C 4 ;g 8 )oo' 1 ' ' 

The K-matrix satisfies the reflection equation, the boundary unitarity relation and the 
boundary crossing symmetry, 

^ 2 (C 2 )^ 2 i(CiC 2 )^i(Ci) J Ri2(Ci/C 2 ) = j R 2 i(Ci/C 2 )^i(Ci) j Ri 2 (CiC 2 )^ 2 (C 2 ), (2.23) 
K{QK{C l ) = 1, (2-24) 

K ?A0 = E<^(C 2 )< Hf 1 C 1 )- (2-25) 

The normalization factor «(£) is determined by the boundary crossing symmetry (|2.25| ). 
The above six relations for the R-matrix and the K-matrix play important roles in solving 
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the model exactly. From these explicit representations those derivatives are 



Q 



d_ 



C=i 



C=i 



1 — r 



a*a* +1 + rrfrrj. , + A<tJ<tJ +1 ) + const., (2.26) 
a z + const., (2.27) 



1 - q 2 v ' ' 
1 + r 



where a x , a y , a z are the Pauli matrices and the anisotropic parameter A = (q + q~ 1 )/2. 
Pi j is the permutation operator: Vi®Vj — > Vj®Vi. In the paper we work in the antiferro- 
magnetic region A < — 1 or — 1 < q < 0. Taking notice of the relation Rj ,j+i(l) = Pjj+i 
and Rjj + i(QPj y i = Pj t iRij + x(Q, we get the derivative of T(£i, . . . , (n) as follows, 



d 

^•(Ci,...,Cat, 



Ci=-=Cw=i 
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^--Ri+i,i(o)-Pi+i,i 
i=i °v 
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+ ^ -p^rRi,i+l{Cj)P,i+l 
=1 i=1 °"->j 



0=i 



N-l 



d 



i=j 
Aq 



Hxxz + const.. 







7V-1 



d 



o=i i=j 



1 -g 2 

Here Hxxz is the Hamiltonian of the XXZ model with boundary magnetic fields 

iV-l 

#xxz = -- E(«+i + °Mu + A « + i) + ^ + ^ 

i=l 

and the boundary magnet fields hi, are related with parameters f, r, 

1 — q 2 1 +f 1 — g 2 1 + r 



0=1 
(2.28) 



IV' 



(2.29) 



4q 1 — r ' 



4q 1 — r 



(2.30) 



We thus have shown that the operator ( |2.13 ) generates the transfer matrix for the finite 
XXZ spin chain. 



To get this transfer matrix we modify the relation ( |2.6|) into 



(2.31) 



We call the states \W) and (V\ the boundary states. Then the boundary q-KZ equation 
is written as 

7;(Ci,...,0v)|Ci,---,Civ> = |Ci,---,0v>, (2.32) 
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with the state . . . , (n) being defined by 



|d,...,Gv> = (vmCi)®---^HCN)\w) 

= (V\<t> ei (t 1 )---$ eN ((N)\W)(v ei ®---®v eN ). (2.33) 



ei,...,eN 



Note that the state • • • , (n) is in an iV-fold tensor space V\ ® • • • ® V^. We regard 
this state as the eigenstate of the transfer matrix of the XXZ model. 

We point out a similarity between such construction by the boundary q-KZ and the 
matrix product ansatz in the one- dimensional react ion- diffusion processes with bound- 
aries flnj. In the matrix product ansatz the eigenstate is expressed by a product of the 
Zamolodchikov-Faddeev algebra, which is determined by the "Hamiltonian" describing 
the time-development of the model [ 2(| . In the viewpoint of the matrix product ansatz, 



the VO is a kind of the Zamolodchikov-Faddeev algebra, since it obeys the same com- 
mutation relation (|2.8|) . However the original matrix product ansatz is not applicable to 
the XXZ model, because the structure of the energy levels may be different. Details are 
discussed in the next section. Our method to construct the eigenstates is restricted to 
the integrable models whose correlation functions satisfy g-difference equations like the 
q-KZ equation. 



3 Energy Levels and Turning Symmetry 

Strictly speaking, there exist no boundary states satisfying the defining relations ( |2.5| ) 
and ( j2.31| ), which include some extra factors A w ((; r) and A v ((; r). These factors deter- 



mine the vacuum energy of the model. In this section we slightly modify the boundary 
states from a mathematical standpoint and also construct the excited states. 

The quantum affme algebra U q (sl2) has two types of the VOs. The VO used in 
the previous section $(C) is the type I, which is defined by the intertwiner for a map 
. v(Ai) -> V(A a _i) ® V (i = 0, 1). Here V(Ai) is the highest weight represen- 
tations of U q {sl2) with the level 1. This mathematical definition leads to 

|W0 (i) e V(Ai), {i) (V\eV*(A t ), (3.1) 

with the dual highest weight representations V*(Aj). Using scalar factors A^(£;r) and 
Ay\(; r) we generalize the relations (|2.5| ) and ( |2.31|) into 

XX(C)*?^ 4) (C)|W0 W = A$((;r)^\C- l Wf\ (3.2) 

e' 

^W(V|$f-)( C -i)^,(C) = A^iar^iVl^^iO- (3.3) 

These scalar factors are given by (see Section |5.2|) 

A^(C;r) = A«(C;r) = 1, A$({; r) = A^(C; r) = A(C; r), (3.4) 
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where 

A(C ' r) -C 2 e 94 K- 2 )e g4 (W)- {6 - b) 

Using the boundary states |W)^ and ^(V|, we redefine the eigenstates of the transfer 
matrix by 

Ki,...,0v) w = ^\v\& i > 1 - i \( 1 )®---®& 1 - i ^(( N )\w)V 

= {i HV\^ l - l \Ci)---®t- l ' t \(N)\Wf\v ei ®---®v eN ). (3.6) 

ei,...,ejv 

The labels (1 — decorating the VOs are omitted if there is no danger of confusion. 

Suppose that the number of sites N is even. We turn the eigenstates upside down, 
which means that the spins and the their orders are reversed. With this operation we 
have 

W^ln^) • ■ ■ ^(CiWir- 1 ))®^ ® • • • ® v_ ei ) 

~ ®(V(r- 1 )MCi)®---®*(CN)\W(f- 1 ))®. (3.7) 



€i,...,e N 



Here we have denoted the parameter r dependence explicitly. Pay attention to the 
defining relations of the boundary magnetic fields ( [2.30|) : the parameter r is inverted 
when the direction of the magnetic field is reversed. If the eigenstates are equivalent to 
the states turned upside down, the relation f = r~ l is derived. We call such geometric 
symmetry the turning symmetry. Since the spontaneous magnetizations are generated by 
the states satisfying the turning symmetry, we hereafter consider only the states whose 
number of sites is even and which satisfy the turning symmetry. Then the eigenvalue 
A (Cii r ) °f ^e transfer matrix Tj(Ci, ■ ■ ■ , Ov) with the eigenstate . . . , (n)^ are 

AC°)(C; r) = A(C; r" 1 ) = A(C; r)~\ A«(C; r) = A(C; r). (3.8) 

Noting that Tj(l, . . . , 1) = A«(l; r) = 1, differentiating T^i, ■ ■ ■ , (n) and A«(C; r) at 
Ci = • • • = Cn = 1, we have the following energy spectrum, 

£r XX z|l,...,l> w = (-l) 1 - < e(r)|l,...,l> w , (3.9) 

2Kk' ixK' 

■ sinh — — sa{2K'a) (-1 < r < 0) 



e r 



2# 1 , , (3-10) 

Here k' is the complementary modulus of the elliptic function and a common constant 
factor is omitted, r and q are parameterized as 

I {q 2 T (0<r<l) ' g ~ 6 • [6 - LL) 
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From this energy spectrum the vacuum state is determined as follows, 



■ 1 |Ci, - - - , Cat> (1) (|r| > 1). 



To assure that the state fl3.12| ) really defines the vacuum, we show that any excited states 
can be constructed by the type II VO ^(C) (/•* — ±) i n the same way as the infinite 
XXZ spin chain |I|, H]. We introduce the M-particles excited states by 



|Ci) • • • > Cat; £i, ■ • ■ , £m) 



w (vi$(Ci) ® • • • ® ^(Civ)*;^) • • • (3.i3) 



and confirm that this state corresponds to the excited state. Using the communication 
relation between the type I and type II VOs, 

*e(0M0 = t(C/0*m(0<MO, (3-14) 

t(C) = 7 2^ ' = (z;P)°o(pz _1 ;p)oo(p;p)oo, (3.15) 
we obtain the eigenvalues, 

^j'(Ci) • • • j Cat) I Ci j • • • j Cat; 6i, • • • , €m) 

M 

= II ^/&MC;&)iCi, • • • , Ov; 6, • • • , Zm) (i) . (3.16) 

i=l 

Then the energy levels of the excited states are written as 

H xxz \l,...,l^ 1 ,...X M )^ = ^\r;^...X M )\l,...,l^ 1 ,...X M f\ (3.17) 

M 

We have parameterized £ = —ie ld . Such a structure of the energy levels is similar to the 



+ V sinh-— dn OA. 3.18 

i=l V 7 



quasi-particles structure pointed out by Faddeev and Takhtajan in the XXX model [£L7 
We conclude that our method gives the vacuum and the excited states. 



4 Thermodynamic Limit 

We have constructed the eigenstates of the XXZ model satisfying the turning symmetry 
in terms of the VOs. With the aid of their free field realizations one can compute 
their 2 Ar -dimensional vector expressions. However they contain the complex integrals 
over N variables, because a component of the VO $+(C) is represented in the integral 
form (see the definition (|5.8| )). To avoid this difficulty, we notice that, in the bulk 
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part, the boundary q-KZ equation requests only the commutation relation of the VOs 
( p.8| ). This fact implies that the vacuum state of the model can be obtained without 
the representation theoretical meanings of the VOs. In this section we decompose the 
vacuum into a huge product of the R-matrices and find its asymptotic form applying the 
Baxter's formula on the corner transfer matrix |T], [Tj| . 



4.1 Half-Infinite Chain Limit 

To calculate magnetizations we begin with the dual eigenstates of the transfer matrix. 
We write them as 

«<0v,...,Cil = w (^|$*(Civ)®---®^(Ci)l^) w 

= Yl (i) (^i $ : JV (a)®---®c 1 (Ci)in {,) (^ ® ••■««*). (4.i) 

ei,...,£iV 

where the boundary states and |V*)M are defined by the relations, 

^(W^AC 1 )^) = A$(C;r)W<W|*:(C), (4.2) 

e' 

E^'(cR*(om (i) = A{?(c;r- i )*:(r i )in w - (4-3) 

e' 

These dual eigenstates form an A-fold dual tensor space <g> • ■ ■ ® V{ and satisfy 

« (Ci, ■ ■ • , CvPKCi, ■ ■ ■ , Car) = A (i) (0; r) <*><&, . • . , Cn\- (4.4) 

Let us explain why the dual states are given by the equation Q4.1). Consider the half- 
infinite chain limit. One of the boundaries then vanishes and the eigenstates are redefined 
by 

|...,Gv,...,Ci> W = ■••®<l>(a)®---®$(Ci)l^) (i) , (4.5) 

«(G,...,Gv,...| = (i) (W*|$*(Ci)®---®r(Civ)®---. (4.6) 

As an example we introduce an operator O on site m whose matrix elements are defined 
by 0.v € = ^2 e , O e e v e i and calculate the vacuum expectation value of O. Use the unitarity 
relation of the VOs ^ e $*(£)$ e (£) = id Jl] and take care of the order of the A-fold tensor 
product of two-dimensional vector spaces. Then the vacuum expectation is written as 

(£>)M = w (wi$: 1 (Ci)--^L(cjo£^(cj---^ 1 (Ci)iw r > (i) - (4.7) 

This result is the same as that in the half-infinite XXZ spin chain except for the nor- 



malization 13 
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4.2 Baxter's Formula 



We consider the thermodynamic limit of the vacuum. We must deal with the infinite 
product of the VOs. Note that, in the bulk part of the vacuum state, only the commu- 
tation relation Q2.8|) has been used. This relation is obtained from the VO represented 
as the infinite product of the R-matrices. Hence the VOs in our method are replaced 
with a product of the R-matrices: 

• • • ® $(Cb) ® $(Cb) ® = ® (•• • %b)%b)) ® (• • • R(( B )R(( B )) <g> • • • , (4.8) 

where -R(C) = R(C)P and Cb stands for the spectral parameter in the balk part. To 
treat such a huge group of the R-matrices we employ the idea of corner transfer matrix, 
which was originally introduced by Baxter to calculate the correlation functions flnj. It 
is defined by 

A(Cb) = R((b) ® (R(Cb)R{Cb)) ® • • • ® (R(Cb)R(Cb) ■ ■ ■ R(Cb)) ® • • • ■ (4.9) 

Baxter discovered that the corner transfer matrix is drastically simplified in the infinite 
chain limit; 

A((b) = Cb°- (4-10) 

Here D is an operator independent of Cb and has the property £,~ D $(()£ D = ( see 
Chapter 4 in [0]). This asymptotic identify is derived from the Yang-Baxter equation 
( [2. 18] ) and also understood in a framework of the crystal base pl| , |22|| . The identity 



Q4.10Q is referred as the Baxter's formula. We remark that the ground states make 
no contribution to the Baxter's formula. Here the ground states mean the completely 
antiferromagnetic states where + and — alternatively appear at each edges of the R- 
matrices. 

By virtue of the Baxter's formula we can simplify the infinite product of the VOs. 
The R-matrices of the VO at the infinitely far point are in the ground states [p], fLlfl . 
In the same manner the R-matrices belonging to the VOs at sufficiently far sites are 
certainly in the ground states. This reasoning implies that the upper triangle sector of 
a huge group is in the ground states. Notice that the R-matrices in the ground states 
are omitted. Then the Baxter's formula gives 

• --HCb) ® $(Ci?)$e(O|W0 (l) = ( B D M(W) {1) . (4.11) 



The VO at a boundary remains because the defining relation (|3.2| ) requests the explicit 
expressions of the VOs. 

Next we consider the product of the VOs near another boundary. Recall that the 
vacuum state satisfies the turning symmetry. Such state is realized by the product of the 
R-matrices and the 90°-rotated matrices. Therefore the turning symmetry tells that the 
product near W(V| is related to that near |W)W through the crossing symmetry ( 2.20|) . 
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Thus, the vacuum state in the infinite chain limit is simplified as follows, 

| C) C B }« = W(K|---®$(-g- 1 C B )®$(-g- 1 Ci ? )$-e(-g- 1 C) 
= (l) (v\(-q- 1 C B 1 r D ®-.(-q- 1 C)CB D MC)\w)V 

= ^(V\$^((( B )(-q) D MO\W)V. (4.12) 

By the definition the spectral parameter ( B of the rotated VOs remains unchanged [l|]. 
Similarly we get the dual state of the vacuum in the infinite chain limit, 

{i) (CXb\ = il HW*m(( B )(-q) D $*_ e (C)\V*) il) . (4.13) 

This infinite chain limit is regarded as the thermodynamic limit because boundaries on 
both sides survive. We thus have simplified the expressions of the vacuum states in the 
thermodynamic limit. 

5 Spontaneous Magnetization in the Finite XXZ 
Spin Chain 

5.1 Free Field Representations of the Vertex Operators 

To calculate the spontaneous magnetization at a boundary, we prepare the free field 
representations of the VOs and realize the boundary states. 

The highest weight representations of U q {sl2) with level 1 are realized as the Fock 
spaces, which consist of the bosonic operators a k {k G Z/{0}) and the zero-modes Aj, a 
(2 = 0,1), 

V(Ai) = C[a_ 1 ,a_ 2 ,...]|0)®(© neZ e A ' +na ), (5.1) 
V*(A,) = (0|C[ai,o 2 ...](8)(© n6Z e- A< - nai ). (5.2) 

Here the bosonic operators satisfy the following commutation relation, 

[a k ,ai\ = — - — d fc +i,o, N = r- ( 5 - 3 ) 

k q — q 

The action of the zero-mode is defined by 

e \ e ? = z d .e? = z^V, (5.4) 

for V(Ai) and 

e?>. = e?.z 9 = e?z^ 9 \ (5.5) 

for V*(Aj). The zero-mode Aj is related with Ai = Ao + a/2. The operator d is 
determined by the commutation relations [d, a] = 2 and [d, A ] = 0. The zero-modes 
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Aj and a are regarded as the fundamental weights and one of the simple roots of SI2, 
respectively. We write 

|0)« = |0)®e A % ®(0\ = (0\®e- Ai . (5.6) 
The VOs of type I are bosonized as follows, 

<I.' '"'(c) = e p «V« 2 )®e«/ 2 (-<7 3 C 2 )^C l , (5-7) 
_ I dw (!-? 2 )^C . -(!-*,<)/ 



/ c 27TV-1 g(w - <rC )(w - Q C ) 



IH = e fi(w) e 5(w) ® e - a uT a . (5.9) 

The integration contour C encircles around in such a way that \q 4 ( 2 \ < \w\ < \q 2 ( 2 \- 
The normal-ordering product : : plays the same role as that in the conventional bosonic 
operators. The "primary fields" are 



00 00 



5.2 Boundary States 

We find the boundary states using the free field representations of the VOs. First 
we realize two boundary states and W(W*|. The commutation relations for 

[A,[A,B]]eC, 

e A e B = e l[A,[A,B]] e [A,B] e B e A^ & B g A = ^A^B ^[B,A]^[[B,A],A] ^ (g ^ 

and the formula 

are useful. From the relation ( |3.2| ) with e = — and the relation ( |4.2[ ) with e = +, we 
derive the boundary states |VT}^' ) and as 

(1 00 L 6k 00 \ 

-^MM a -' + ^ / ' ,a - t |0>< "'' (5 ' 14) 
k — 1 k — 1 / 



-2k 



00 



1 \ ^ kq~ 



(in = «)<0|exp + . (5.15) 
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Here the functions fj? and are given by 

/« = ^ H,(l - q k ) + (-iJ^gd-aOV 1 -^*) , (5.16) 

= ^(^(l-rt + C-l) 1 ^- 30 *^ 1 -*)*), (5.17) 

with 6k = (1 + (— l) fc )/2. For the relations with other components we must verify these 
realizations. As an example we consider the relation ( |3.2| ) with e = + in the case i — 0. 
The realization ( |5.14|) reduces the relation (|3.2| ) to the following integral equation, 

w-'dw (l-r( 2 )(w 2 -q 6 )(w-rq 4 ) (0) = r -i^i n \(o) 

2 7rv /3T( w -g2 C 2)(^- (? 4 C 2 )(^- (? 4 C- 1 ) 1/ ^ ^ / • 

(5.18) 

In the left hand side we change integration variable as w — > q e w 1 and take the average of 
the original and the transformed terms. We then get the symmetrical integral expression 
for the spectral parameter ( and 

Other boundary states are found from the following identities, 

= (^vTK-g)- , \y*) {l) = (-q)~ D \W)®. (5.19) 

They are proved by making use of the equation Kf((;r) = A((;r)KZ^ ((;r~ l ). For 
instance, putting the first identity of ( |5.19| ) into the relation (|4.2|), we have 

e' 

= A((;r) ^W^K-g^h-^/Q^r- 1 ) 

e' 

= A(C;r)A«(C;r- 1 ) w (V|$_ e (C)(-g) D 

= Ag(C;r)«(W*|$:(C). (5.20) 



The defining relations ( |3.2| ) and ( [4.2| ) are reproduced by the relations ( |4.3j ) and ( |3.3| ) 
through the identities ( |5.19| ) respectively. 

We explain briefly why factors A^-((;r) and Ay\(;r) appear in the defining rela- 



tion of the boundary states. In section |4.1| we have formulated the dual eigenstates by 
comparing the correlation function in our method with the already known result. This 
argument is valid only if the states and W(W*| are the completely same as those in 

the half-infinite XXZ spin chain. Thus, factors A^(£; r) are attached to these states (see 
|T3|). The factors in all defining relations are determined through the relations (|5.19|) . 



5.3 Spontaneous Magnetization 

We are in a position to calculate the spontaneous magnetization in the finite XXZ spin 
chain. By the Baxter's formula ( |4.10| ) and the identities ( |5.19|) , the vacuum states are 
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written as 

(c±) (i) = «(V|$_ e (C-)(-g) D $ £ (C + )l^)« 

= «(iy*|$_ e (-g- 1 C-)$ e (C + )|^)«, (5.21) 

«( C± | = «(ty*|$:(c_)(- g ) D $*_ e (c + )|y*)« 

= w (^*|$_ e (-g- 1 C-)<l'e(C+)|W / ) W - (5.22) 
Therefore only the following quantity is of our main concern, 

PiH&r) = W(W|$_ e (-g- 1 C-)$ e (C+)|W0 W - (5-23) 
Noting the order of the sites in the dual state we have the magnetization at a boundary, 

,-.),«,.) (if(l;r))*-(^>(l;r))3 

By the free field representations of the VOs we get the following bosonic expression of 

Pe w tet;r), 

±,r e (g 6 ;g 4 )oo Jc« 2^^ (w - *_)(«; - g 2 z £ )(w - g 4 ^+) 

The integration contours and C_ encircle clockwise around in such a way that 
|g 4 z±| < \w\ < \q 2 z±\ and \q 2 z±\ < \w\ < \z±\, respectively. Y®(z±,w;r) is the inner 
product of normal-ordered bosonic operators given by 

Y®(z±,w;r) = ^(W*\exp J exp ( -f^a k F k ) \W)®, (5.26) 



v fc=l / \ k=l 



where the functions Fu and Gu are 



ft = V ( ^ +?2 ^ l) "V"' (5 - 27) 

7fc/2 fc/2 

G * = w (4+, " V) "V m '- (5 ' 28) 



Since both boundary states ( |5.14|) and (|5.15|) contain the quadratic terms of bosonic 
operators on the exponential function, the commutation relations among bosons are not 
useful. First we insert twice the completeness relation, 

idvw = / m 2 ^V (-E|^l^| 2 ) \0 {m (ii (5-29) 

between the boundary states and the exponential operators. Here the integration is 
taken over the entire complex plane and the measure T>\^\ 2 means the infinite product 
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HfcLi^l^fcl 2 except for the normalization. The coherent states and are given 
by 

l«> ( " = ^(gj^&^lOP, (5.30) 

They clearly satisfy the definition of the coherent states; 

a fc |0 W =&IO ( V°<fla-fc = W <el6fc. (5-32) 
Next we use the formula for the Gauss integration, 

p / oo i oo 

'^M!M W y " (5 .33) 



cxp 



fc=i 



For Y^(z±,w;r) is an invertible 4x4 matrix and is a four-component vector. 
Then the inner product Y ( - t \z±,w;r) is reduced to 

V i — i 9 - 



'{9k ~Q Jk )fik + {Jk - Q 9k J G . 

1 -2k Ai)2 Ai)Ai) , * 6fc„(i)2 



+§?-*/r-/w+§^nj- (5 - 34) 

Taking out terms depending on the parameter e, z and tu, we have the following integral 
expression, 

p(i) (z . r) =I edw !-^^" 2 ®Aq 2 w 2 ) (535) 

£ 1 ' ' J c ^27r^l-wrq- 2 Q q2 (wz)e q 2(wz- 1 ) { } 

Here we have put z = z + = Z-, since we have no poles that degenerate with this identi- 
fication. The integration contour picks out the points q 2n z, q 2n z~ x for n — 2, 3, . . . 
and the contour the same points for n — 1,2,.... The contour C± encircles around 
not only the points in C±^ but the point q 2 r~ x . With the aid of the quasi-periodicity 
property Q p (pz) = —z^ 1 Q p {z) the integral expression (|5.35| ) is easily calculated. We 
thus arrive at the following results: 

(5.36) 



+ V- 1 ' 

? M (l;r) = -(l+Pj 0) (l;r)), P e (1) (l; r) = -P e (0) (l; r" 1 ). (5.37) 
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Here we have normalized the result such that the sum of P + and P_ equals —1. With 
P«(l;r) the magnetization at a boundary is given by 

(v*)W(r) ~ l + 2P(0(l;r) 

(a) [r) ~ l + 2P(*)(l;r) + 2(PW(l;r))2- 

We see that the magnetization survives even if the boundary magnetic fields vanish. 
Substituting —1 into the parameter r (see Q2.30D ), we obtain the spontaneous magneti- 
zation in the finite XXZ spin chain as follows, 

The spontaneous magnetization in the half-infinite XXZ spin chain is known to be 
-(9W)^c/(-<?W)L @. Denoting it by (a 2 )^ we have 

(a z > = 2<o*>oc " 2(OL + 2<OL - • • • • (5-40) 

The relation ( 5.40| ) shows that the spontaneous magnetization in the finite chain is 
roughly twice as large as that in the half-infinite chain when q is close to —1 or when 
the anisotropy is small. We conclude that, in the XXZ model, the anisotropy generates 
the spontaneous magnetization and the finiteness of the model enhances it. 



6 Summary and Discussion 

We have studied the finite XXZ spin chain with boundaries. The vacuum state is derived 
from the boundary q-KZ equation and is composed by the type I VOs. Apart from the 
representation theory we have interpreted a large product of the VOs as a huge group of 
the R-matrices and have simplified the vacuum state in the thermodynamic limit using 
the Baxter's formula. Simple expressions of the states enable us to calculate the spon- 
taneous magnetization in the finite XXZ spin chain. The spontaneous magnetizations 
in the finite chain and in the half-infinite chain are related by ( |5.40| ). It has been made 



clear that the vertex operator approach is valid for the analysis of the finite XXZ spin 
chain. 

We have constructed the excited states by including the type II VOs in the vacuum 
and have shown that the energy levels are similar to the quasi-particle structure in the 
XXX model ||17|| . Such structure is known to appear in the infinite chain limit. We 
bypass arguments whether there exist these energy levels in a finite chain, since we 
attach importance to the thermodynamic limit. 

In the paper we have considered only the diagonal K-matrix. It is interesting to solve 
the same problem when the K-matrix has non-diagonal elements. The non-diagonal K- 
matrix generates not only the boundary magnet fields but the boundary flows of spinons 
f23fl. Since the sum of spins is not preserved we might expect a novel phenomenon. To 
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treat the non-diagonal K-matrix one needs the complex integrals, because each defining 
relation of the boundary states is expressed by an equation such that two VOs are 
combined into one. Such equations can not be realized in terms of only bosonic operators. 

We have regarded the infinite product of the VOs as the thermodynamic limit since 
both boundaries remain. The thermodynamic limit should be taken by fixing a ratio of 
two infinite parameters. In our method we may choose the ratio between the number 
of the type II VOs and the number of the sites. In the paper, however, we have not 
given attention to the fixed parameter, because only the vacuum state is considered. In 
return this choice implies that one must treat the infinite number of the type II VOs 
when one compute the thermodynamic limit of the excited states. It is difficult to find 
the behavior of the infinite product of the type II VOs by their free field representations. 
A way to solve the difficulty is to find the recursion relations among the excited states. 

Korepin proposed a method to calculate the correlation functions in the integrable 
systems |24|. There, the recursion relations are derived from the algebraic Bethe ansatz. 



There may exist a similarity between the recursion relations by our method and those 
by the algebraic Bethe ansatz. Referring to this similarity one can realize the Korepin's 
method by the vertex operator approach. This expectation is supported by a fact that 
some combinations of the type I VOs and the type II VOs correspond to the L-operator 
in the Yang Baxter equation p5| . 
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